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Abstract: We consider a Higgs boson coupled to gluons via the five-dimensional effective 
operator HtTG ^^G^^ . We treat H as the real part of a complex field 4> that couples to 
the selfdual gluon field strengths and compute the one-loop corrections to the (j)-MSM 
amplitudes involving 0, two negative helicity gluons and an arbitrary number of positive 
helicity gluons. Our results generalise earlier work where the two negative helicity gluons 
were constrained to be colour adjacent. We use four-dimensional unitarity to construct the 
cut-containing contributions and the recently developed recursion relations to obtain the 
rational contribution for an arbitrary number of external gluons. We solve the recursion 
relations and give explicit results for up to four external gluons. These amplitudes are 
relevant for Higgs plus jet production via gluon fusion in the limit where the top quark 
mass is large compared to all other scales in the problem. 
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1. Introduction 

The startup of the LHC anticipated for the autumn of 2008 heralds the arrival of a new 
arena for the exploration of particle physics. The large centre of mass energy is expected to 
produce complex multiparticle final states both as decay products of putative new physics 
Beyond the Standard Model and through the Standard Model itself. Extracting the signals 
of new phenomena and discriminating between different models of new physics is only 
possible if the predictions for the Standard Model, and its prominent extensions, have 
sufficient accuracy. The precision which can be achieved using calculations at leading 
order in perturbation theory is, in most cases, not sufficient for detailed studies of signals 
and especially backgrounds at the LHC. In many cases, the calculation of multi-particle 
final states at next to leading order (NLO) will be essential to the successful interpretation 
of the data. Over the past few years vast leaps in our understanding of the structure of one- 
loop amplitudes in gauge theories has lead to the widespread belief that soon predictions 
for many multi-jet final states will soon become available. 

The use of four-dimensional on-shell techniques, originally pioneered by Bern ct al [1, 2] 
in the mid-90's has lead to a vast reduction in the complexity of one-loop calculations. The 
use of gauge-invariant physical amplitudes (at tree level) as building blocks means that 
simplifications due to the large cancellation of Feynman diagrams occur in the preliminary 
stages of the calculation, rather than the latter. The unitarity method sews together four- 
dimensional tree-level amplitudes and, using unitarity to reconstruct the (poly) logarithmic 
cut constructible part of the amplitude, successfully reproduces the coefficients of the cut- 
constructible pieces of a one-loop amplitude. This has extensive uses in supersymmetric 
Yang-Mills theories, which are cut-constructible i.e. the whole amplitude can be recon- 
structed from knowledge of its discontinuities. 

The more modern applications of unitarity were kick-started by the discovery of the 
MHV rules by Witten and collaborators in 2004 [3]. The realisation that MHV tree am- 
plitudes could be promoted to vertices which could be used to create amplitudes with any 
number of negative helicity gluons sparked a revolution in the field of on-shell QCD. In a 
series of remarkable papers, Brandhuber, Spence and Travaglini (BST) [4] showed how the 
MHV rules can be used at one-loop for the calculation of n-point gluonic MHV amplitudes. 
Around the same time, the quadruple cut [5] using complex momenta was introduced to re- 
duce the determination of the coefficients of box integrals to simple algebraic manipulation 
of four tree level amplitudes. Double and triple unitarity cuts have led to direct techniques 
for extracting triangle and bubble integral coefficients analytically [6-9]. In cases where 
fewer than four denominators axe cut, the loop momentum is not frozen, so the explicit 
integration over the phase space is still required. In the BBCFM-approach [6-8], double 
or triple cut phase-space integration has been reduced to extraction of residues in spinor 
variables, and, in the case of a triple cut, residues in a Feynman parameter. This method 
has been recently used for the evaluation of the complete six-photon amplitudes [10, 11]. 

Despite its success, the four-dimensional unitarity method does not give the complete 
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result for non-supersymmetric theories such as QCD, since there are missing rational func- 
tions which are cut-free and as result do not possess discontinuities in physical channels. 
The missing rational parts have only simple poles and are therefore tree-like. Since the 
rational pieces of one-loop amplitudes are tree like in their discontinunity structure they 
can be calculated using a straightforward generalization of the tree level recursion relations. 
One can then use the tree-level on-shell recursion relations [12, 13] to compute the rational 
pieces of one-loop amplitudes recursively. The ability to calculate the rational pieces of 
amplitudes independently of the cut-constructible terms lead to the development of the 
unitarity bootstrap approach [14-20]. Recently, an automated package BlackHat has been 
developed to compute these rational terms for pure QCD amplitudes [21]. 

Another approach is to extend use of unitarity to D = 4 — 2e dimensions [8, 22-27] 
and to take the cut particles into D = 4 — 2e dimensions. This approach has the great 
advantage of calculating both the cut containing and the elusive rational terms at once, 
but care must be taken with application of the four-dimensional spinor helicity formalism 
in D dimensions. 

It has also been observed that the rational parts are related to the ultraviolet behaviour 
of the amplitude, and can be directly obtained from the traditional Feynman diagram 
approach [10, 28-30]. In a very interesting work, Ossola, Papadopoulos and Pittau [31] have 
applied the unitarity ideas directly to the integrand of the Feynman amplitude, without 
necessarily appealing to the simplified forms of the cut diagrams. They find algebraic 
identities which can be automatically solved to give the coefficients of the master integrals 
as well as the rational part. This approach is being further developed [32-38] with a view to 
providing automated computations of both cut-constructible and rational parts of one-loop 
scattering amplitudes. A summary of the current state of the art is given in Ref. [39]. 

In this paper, we exploit the unitarity bootstrap approach [14-19] which meshes to- 
gether the calculation of the cut-constructible parts of an amplitude (via generalised uni- 
tarity, one-loop MHV rules etc.) with the ability of the BCFW recursion relations to 
calculate the rational pieces. As a result of the splitting the total amplitude is given by 
the combination 



Here the Cn are the purely cut-constructible pieces which arise from box, triangle 
and bubble (and in massive theories tadpole) loop integrals, the functions in Cn are those 
which contain discontinuities, in general poly-logarithims (and associated tt^ terms). C„ 
may contain unphysical singularities which are produced by tensor loop integrals and must 
be cancelled by rational contributions. To make this cancellation explicit, we add the 
cut-completion terms CRn, so that the "full" cut-constructible pieces are defined as, 



(1.1) 



n 



(1.2) 
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These additional rational terms would be double counted if we naively calculated the 
rational terms with the BCFW recursion relations, so we redefine the rational pieces as 



Rn — Rn — CRn- (1.3) 

The rational part now contains only simple poles, and can, in principle, be constructed 
recursively using the multiparticle factorisation properties of amplitudes. We label this 
direct recursive term by . By construction, the recursive approach generates the com- 
plete residues of physical poles. However, the cut-completion term CRn may also produce 
a contribution at the residue of the physical poles, and may lead to double counting. These 
potential unwanted contributions are removed by the overlap terms, 0„. 

To generate the recursive contribution, one generally shifts two of the external momenta 
by an amount proportional to z. Complex analysis [13] then generates the correct amplitude 
provided that 

An{z) ^0 as z ^ oo. (1.4) 



For a generic tree- level process it is frequently possible to shift two momenta such that ( lA ) 
is obeyed. Similarly, for one-loop processes, one can often make a similar shift. However, 
because the choice of CRn is not unique, the shift may introduce a "spurious" large z 
behaviour in CRn, labelled by Inf Cii„, which should be explicitly removed [18, 19]. The 
rational part (provided that An{z) — > as z ^ co) is given by, 

Rn = R^ + On-lniCRn, (1.5) 

while the physical one-loop amplitude is given by [18, 19], 

=Cn + CRn + R^ + On- Inf CRn- (1.6) 



In this paper, we focus on the 0-MHV amplitudes involving (p, two negative helicity 
gluons and an arbitrary number of positive helicity gluons. Our results generalise earlier 
work [20] where the two negative helicity gluons were constrained to be colour adjacent. 
The paper proceeds as follows. In section |2|, we give a brief overview of the Higgs couples 
to gluons, and how this is related to (/(-amplitudes. Section ^ reviews the four-dimensional 
unitarity methods for constructing the cut-containing contribution C„. There are many 
similarities with the pure-gluon case, and we develop the derivation of the cut-constructible 



parts of pure-glue MHV amplitudes and (/>-MHV amplitudes in sections 3.S and 3.4. Sec 



tion ^ deals with computation of the three separate rational pieces, the cut-reconstructible 
part CRn, the on-shell recursive part Rn and the overlap term On- As an example, we 
derive the four-point amplitudes ■*((/), 1", 2+, 3^, 4+) and yl^^^ (//, 1^, 2+, 3~, 4+) in sec- 
tion |5|, while section |6| describes the checks we have performed on our result. Finally, in 
section |7[ we present our conclusions. Two appendices detailing the explicit construction 
of the cut-completion terms and the forms of the one- loop basis functions are enclosed. 
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2. The Higgs Model 



The coupling of the Higgs to gluons in the Standard Model is produced via a fermion 
loop. Since the Yukawa coupling depends on the mass of the fermion, the interaction is 
dominated by the top quark loop. For large mt this can be integrated out, leading to an 
effective interaction, 

= ^iftrGp^G^^ (2.1) 

This approximation works very well when the kinematic scales involved are smaller than 
twice the top quark mass [40-42]. For the interesting pp ^ H plus two jet process, the 
approximation is valid when mn, PTj < nT-t [43]. The strength of the interaction C has 
been calculated through to order ©(af) in the standard model [44]. To order 0{al) [45], 
this is 

C = ^(l+'-^^ + ...] (2.2) 
bnv \ 4 vr 

The MHV-structure of Higgs-plus-gluons is best understood [46] by defining the Higgs 
to be the real part of a complex scalar (f> = ^{H + lA) so that 



(2.3) 



where the purely selfdual (SD) and purely anti-selfdual gluon field strength tensors are 
given as 

Gg'^ = ]^{G^"' + G'XsD = liG^" - (2-4) 



with 



'G"" = ^e'^^^^Gp^. (2.5) 



Because of selfduality, the amplitudes for (j) and (f)^ have a simpler structure than those for 
the Higgs field [46] . The following relations allow for the construction of Higgs amplitudes 
from those involving (j) and (j)^ . 

AtHH,g^\...,g^") = Alr\<P,9^\...,9t)+At\<l>l9^\...,9t), (2.6) 

4"^) ( A, 5^ , . . . , 5> ) = ^ (4"^) (</., 5^ , . . . , 9^ ) - ^^T^ (0^ . . . , )) . (2.7) 

Furthermore parity relates (j) and (j)^ amplitudes, 

{<p\g^\..., 5^ ) = (4"^) (0, 5r'^ , • • • , 9n'" )J- (2.8) 

From now on, we will only consider 0-amplitudes, knowing that all others can be obtained 
using eqs. (|^)-(U). 
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The tree level amplitudes linking a (j) with n gluons can be decomposed into colour 
ordered amplitudes as [47, 48], 

A^^\<j>,{ki,\i,ai}) = iCg^-^ tr(TMi)...rMn))40)(^^^(lAi^__^„A„))_ (2.9) 

a£Sn/Z„ 

Here Sn/Z^ is the group of non-cyclic permutations on n symbols, and j^^ labels the 
momentum pj and helicity Xj of the j^^ gluon, which carries the adjoint representation 
index Oj. The T"' are fundamental representation SU(iVc) color matrices, normalized so 
that Tr(r«r^) = 5"^. The strong coupling constant is = c/VC^tt). 

Tree-level amplitudes with a single quark-antiquark pair can be decomposed into 
colour-ordered amplitudes as follows, 



A^^\<P, fe, Xi,ai}, {Pj,Xj,ij}) (2.10) 
= iCg--^ Yl (^"'^^'^ • • • T''^'"-''h^n M4>, 1\ct{2^\. . . , (n - 1)^""^), n"^) . 

aESn-2 



where S'n-2 is the set of permutations of (n — 2) gluons. Quarks are characterised with 
fundamental colour label ij and helicity Xj for j = 1, n. By current conservation, the quark 
and antiquark helicities are related such that Ai = — A„ = A where A = it^. 

The one-loop amplitudes which are the main subject of this paper follow the same 
colour ordering as the pure QCD amplitudes [1, 49] and can be decomposed as [20, 50, 51], 

[n/2]+l 

A^^\ct>,{ki,Xi,ai})=iCg^ J2 E G'„;e(a)4i)(</.,cT(l^\ • • • ,n^")) (2.11) 

C=l <T&Sn/Sn;c 

where 

Gn;i(i) = Artr(r"i • • • r»") (2.12) 

Gn;c(i) = tr(r"i • • • r"'^-!) tr(r"= • • • t"") ,02. (2.13) 

The sub-leading terms can be computed by summing over various permutations of the 
leading colour amplitudes [1] 

The tree level ^MHV amplitude has the same form as the pure-glue MHV amplitude, 

4O)(0, 1-, 2+, . . . , m-, . . . , n+) = ^^^l^^'^^^y (2.14) 

The only difference between the gluon only and the 0-MHV amplitude being momentum 
conservation, here the sum of all the gluon momenta equals — p</). Since we will encounter 
MHV diagrams in which a fermion circulates in the loop we will also need the amplitudes 
involving a (f) with a quark anti-quark pair [52], 

(1 m)^ (nm) 



ylW((/),l-,2+,...,m-,...,n+) 



(12). ..(nl)' 
(nm)^ (1 m) 
«' (12). ..(nl)' 
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Also as a consequence of the 1-loop nature of the (/)-gluon vertex the following all minus 
amplitude is non-zero at tree-level; 

4 

2-, . . . , n-) = (2.16) 

Amplitudes with fewer (but more than two) negative helicities have been computed 
with Feynman diagrams (up to 4 partons) in Ref. [48] and using MHV rules and on-shell 
recursion relations in Refs. [46, 52]. The MHV amplitude for an arbitrary number of gluons 
but with two adjacent negative helicity gluons was computed in Refs. [20, 53]. 

In this paper we concentrate on the general helicity case for the one-loop (jy-MHV ampli- 
tude. For definiteness, we focus on the specific helicity configuration (1~, . . . , m~ , • • • , n"*"). 
Throughout, we will use the notation, 

Si,j = ijpi + Pi+i H \- Pj-i + Pjf = Pfi^j) 

Sij = 2ipi.p,) = {ij)[ji], (2.17) 
with the exception of section ^ where we use the notation Pabc to represent Pa + Pb + Pc- 

3. The cut-constructible parts 

The calculation of the cut-constructible terms has been performed within both the EST 
approach [4] and the BBCFM approach [6, 7]. Both methods rely on reconstructing the 
amplitude using four-dimensional unitarity with a double cut. Compared to conventional 
methods, one is attempting to compute the (four-dimensional) coefficients of the loop 
integrals as efficiently as possible. The methods differ in how the integration over the 
phase space of the cut particles is carried out. The BST method uses Passarino-Veltman 
techniques to eliminate any remaining tensor integrals, and aims to cast the integrand into 
the form of well-known phase space integrals. It has been shown to work well for MHV 
amplitudes. 

On the other hand, in the BBCFM method, the use of spinor variables yields an al- 
ternative to the Passarino-Veltman reduction of tensor integrals, based on spinor algebraic 
manipulation and integration of complex analytic functions. It has been applied success- 
fully to non-MHV amplitudes. Here, we use both methods as a check of our results. 

3.1 The BST approach 

In the BST approach [4] a generic diagram can be written: 

'^ = ^ I ^^^^'HL,-L2-P)AL{h,-P,-l2)AR{l2,P,-h) (3.1) 
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Figure 1: A generic one-loop MHV diagram or unitarity cut. 



where ^L(R) 3.re the amphtudes for the left (right) vertices and P is the sum of momenta 
incoming to the right hand amphtude. The key step in the evaluation of this expression is 
to re- write the integration measure as an integral over the on-shell degrees of freedom and 
a separate integral over the complex variable z [4]: 

= {z -z')t'+ ,./ hd%5(^\ll)5(^\ll), (3.2) 

where z = zi — Z2 and z' = zi + Z2- The integrand can only depend on z, z' through the 
momentum conserving delta function, 

_L2-P) = {h-l2-P + z'n) = ih-h-P), (3.3) 

where P = P — zrj. This means that the integral over z' can be performed so that, 

^ = ^^^ ^V%<5(+)(Z2)<5(+)(/i)<5(4)(/i-Z2-P)AL(Zl,-P,-/2)A^(/2,P,-Zl) 

= {Aif2Tri J^J dLIPS^''\-h,l2,P)ALili,-P,-l2)AR{l2,P,-h), (3.4) 
where, 

dLIPS(^)(-/i,Z2,P) = -^d%d%d(+\lj)S^+Hq)S^^Hh -I2-P) (3.5) 

The phase space integral is regulated using dimensional regularisation. Tensor integrals 
arising from the product of tree amplitudes can be reduced to scalar integrals either by using 
spinor algebra or standard Passarino-Veltman reduction. The remaining scalar integrals 
have been evaluated previously by van Neerven [54] . 

At this point, one has obtained the discontinuity, or imaginary part, of the amplitude. 
However, by making a change of variables the final integration over the z variable can be 
cast as a dispersion integral 

^ = (3.6) 

that re-constructs the full (cut-constructible part of the) amplitude. 
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3.2 Spinorial Integration 



In the BBCFM approach [6, 7], we make a conventional double cut, so that a generic 
diagram can be written: 

^= (2^/ ^^^^'^i^^-^^-P)Mh-P,-l2)ARih,P,-li), (3.7) 

with ll = ll = 0. 

The double-cut can be written as, 

V = j dLIPS(^) Al{Ii, -P, -12)Ar{12,P, -h), (3.8) 
where the dLIPS^^) can be parametrised in spinorial variables, as follows [3], 

where the delta function eliminates the integration over I2, and the remaining li integration 
variable has been rescaled, li = t I, corresponding to, 

l^i) = Vi\£) , \li] = Vt \e] (3.10) 

with if = £^ = 0. Accordingly, the double-cut can be written as, 

where we indicate only the dependence of the tree-level amplitudes on the integration 
variables. By means of Schouten identities, one can disentangle the dependence on \i) and 
\£\, and express the result of the i-integration (trivialised by the presence of the (5-function) 
as a combination of terms whose general form looks like, 

^= (^E j{^mdl\1i, (3.12) 



with 



{i\ Pi \er+' {i\ P2 \£] 



where Pi and P2 can either be equal to the cut-momentum P, or be a linear combination 
of external vectors; and where the pj's depend solely on one spinor flavour, say \i) (and 
not on and may contain poles in \i) through factors like 1/ (IQ) (with \Q) being a 
massless spinor, either associated to any of the external legs, say \ki), or to the action of a 
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vector on it, like f\ki\). 

The explicit form of the vectors Pi and P2 in eq. (3.13) is determining the nature of the 
double-cut, logarithmic or not, and correspondingly the topology of the diagram which is 
associated to. Let us distinguish among the two possibilities one encounters, in carrying 
on the spinor integration of T^: 



1. Pi = P2 = P (momentum across the cut). In this case, the result is rational, 
hence containing only the cut of the 2-point function with external momentum P (or 
degenerate 3-point functions which can be expressed as combination of 2-point ones). 

2. Pi = P, P2 7^ P, or Pi 7^ P2 7^ P. In this case, the result is logarithmic, hence 
containing the cut of a linear combination of n-point functions with n > 3. 



If Pi = P. = P, 



p\l] 



n+2 



(3.14) 



If, however. Pi = P, P2 7^ P or Pi 7^ P2 7^ P. one proceeds by introducing a Feynman 
parameter, to write Tj as. 



li = {n + l) l\x{l-xr pi{\i)) ^''^^^ 







{t\R\t] 



n+2 



(3.15) 



with 



xf>i + {l-x) 



(3.16) 



The spinorial structure of eq. ( 3.14| ) and eq. ( p. 15 ) is the same. Therefore, we proceed with 
the spinor integration of eq. ( 3. IS]) because it is more general than the eq. ( 3.14 ), because 
of the presence of the Feynman parameter. 

First, the order of the integrations over the spinor variables and over the Feynman 
parameter is exchanged and we perform the integration over the |£]-variable by parts, 
using [6] 



P\£] 



n+2 



(n + 1) {i\ p \£] 



n+l 



P\rj] 



(3.17) 



obtaining, 



(2vr)4 



{I d£)[i di] li 



dx (1 — x) 



^Jiidm 



n+l 



{£\R\£r^^ {l\R\r,] 



(3.18) 



Afterwards, the integration over the |£)-variable is achieved using Cauchy's residues theo- 
rem, in the fashion of the holomorphic anomaly [55-57], by taking the residues at \£) = $\r]\ 
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Figure 2: The MHV diagrams contributing to one-loop gluonic MHV amplitudes 



and at the simple poles of pi, say \t} = \l 



^3h 



1 



(2vr)4 
(27ri) 
(2^ 



J {i di)[i d£] 



dx (1 



\(i?2)n+l ^i^eij^ ^' 



n+1 



R\r^' {i\R\v] 



(3.19) 



To complete the integration of eq. ( 3.1S| ), one has to perform the parametric integration 
which is finally responsible for the appearence of logarithmic terms in the double-cut. Alter- 
natively, the spinorial integration of eq. (3.14) would generate a pure rational contribution. 
We remark that the role of \i) and \i] in the integration could be interchanged. 



3.3 Gluonic amplitudes 



We note that there are many similarities between (/>-amplitudes and pure glue amplitudes, 
and we will exploit this by first rederiving the cut-constructible contribution to pure glue 
MHV amplitudes with the same helicity configuration. 

The graphs contributing to the one-loop gluonic amplitude An\l~ , ■ ■ ■ , m~ , . . . , n+) 
are shown in Fig. ^. There are two distinct types of diagram, labelled (a) and (b). In type 
(a), only gluons circulate in the loop, while in type (b) gluons, fermions (and scalars) may 
circulate. They can be characterised by the following sums 



n—2 n m— 3 m—1 n m—1 n—1 m—1 m m—2 

wEE+EE WEE+ E E + EE- (3-20) 

i=m j=i+2 i=l j=i+2 i=n j=2 i=m+l j=l i=m j=l 

The various contributions have been computed using the MHV rules in Refs. [4, 58, 59]. 
We note that contributions of type (a) associated with a cut in the S(j^;^) j channel have 
an integrand of the form, 

I L R)^j+,),, {h{^ + l))■■■{ji2){i2h){i2U + l))■■■{^h){ili2) 



(12) •••(nl) 

4°) g(i,i + l,j,j + l) (3.21) 



-10- 



where 

(Wi)(i(i + 1)) (fA>0-(i + i)) 

For diagrams of type (b), there are three possible contributions - depending on whether 
gluons, fermions (or for supersymmetric theories scalars) are circulating in the loop. It is 
convenient to consider both (b)-type diagrams in the sq+i) j channel together. Immediately, 
we write down 



{AlAr 
{AlAr) 



(^i^^i + 1)) . . . {j£2){£2h){h{j + !))••• (i^l)(4^2 

fermions 



{lii){U2f{m£2){mhf + {li2){lhf{mh){mhf 



(i+i),* + 1)) . . . {je2){hh)mj + 1)) • • • {iei){til 



21 



/A A N scalars 
[AlAr)(- - 



'U+i),^ ^l^ii + 1)) . . . {j£2){t2h)Mj + 1)) • • • {ih){hh) 



(3.23) 



In each case, the denominator has the same structure as in the (a)-type diagrams and only 
the numerator changes with the particle type. We can exploit the Schouten identity 

(1^2) {mil) - {Ui) {mi2) + (Im) (^1^2) = (3.24) 
to rewrite each of the numerators into a simpler form. 

{U2)^{mh)' + {Ui)'{m£2)^ = {lm)'{h£2)' 

+4{U2){mei){l£i){me2){lmf{h£2f 
+2{U2f{m£if{Uif{m£2f 

(3.25) 

{l£i){l£2f{m£2){m£if + (l£2)(l^i)^(m^i)(m^2)^ = {U2) {ni£i) {l£i) {m£2) {Imf {£i£2f 

+2{l£2f{m£if{l£if{m£2f 

(3.26) 



We see that the first term on the RHS of eq. ( 3.25 ) corresponds to an (a)-type gluonic 
contribution which we label with G, while the third term looks like the scalar contribution 
of eq. ( p.23| ) which we label with S. Similarly, the fermion contribution can be separated 
into a fermionic piece F and a scalar contribution S. We define the three contributions as, 

{lm)^{£i£2)^ 



(ALARf 
(AlAr) 



(^■+1)'* {£i{i + 1)) • • • {j£2){£2£i){£2{j + 1)) • • • {i£i){£i£2) 

= 4°) g(i,i + l,j,j + l) (3.27) 
p _ {l£2){m£i){l£i){m£2){lmf{£i£2f 



(AlAr)^, 



O'+i),* + 1)) . . . {j£2){£2£i){£2{j + !))••• {i£i){£i£2) 

= -4°) .F(i,i + l,j,i + l) (3.28) 
{l£if{l£2f{m£2f{m£if 



O'+i),* (^i(i + !))••• {]£2){£2£i){£2{3 + !))••• {i£i){£i£2) 

= -4°) 5(i,i + l,j,i + l) (3.29) 
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where + + 1) is defined in eq. (|]2|) and, 

. . , _ {^{^ + l)){j{j + l)){U^)^m£,)^le2f{m£,f 
^ ' '-^'-^ ^ (lm)4(£i^2)2(i£i)(^i(^ + l)){ji2){i2U + 1)) ^ ^ 

Restoring the particle multiphcities in supersymmetric theories, we see that for = 4 
SYM with four fermions and six scalars (in the adjoint representation), only the "gluonic" 
part remains 

iALAR)f;:^l^ - 4 (A^A^j;— + 6 (^LAi?)( -S):. = (^l^k)8+i), . (3.32) 

On the other hand, for QCD with Np fermion flavours in the fundamental representation, 
the contribution from this graph is, 

(3.33) 



The functions X for X = G,F, S represent contributions to the cut amplitude. Per- 
forming the phase space and dispersion integrals generates the "cut-constructible" contri- 
bution to the full amplitude. We define, 

X{i,i + l,j,j + l) = I ^ I d''LIPS{-h,l2,P) + + (3.34) 

Explicit expressions for X(i,i + + 1) are written down in Appendix The one- 

loop gluonic MHV amplitude is thus obtained by summing combinations of the "cut- 
constructible" contributions according to eq. (3.20). As a result the one-loop gluonic MHV 
amplitude is given by, 

C„;i(l-,2+,...,m-,...,n+) 

= cMf (^^^i(m,n) - 4(1 - ^^A^^^(m,n) - 2^1 - ^) Af^^(m,n)) (3.35) 

where 

An-i{m,n) = -- ^F4™(si,i+2;si,i+i,Sj+i,i+2) ~ J ^f^'^^Sij, Si+ij-i; Si+ij, Sij-i). 

(3.36) 



i=l i=l j=i+3 



The terms associated with the fermion loop have the following form: 

n m—1 



i=m+l j=2 
m—1 n 



2^ 2^ ^2 -^lm-'l(^(i+l,i)'^(i,j)) 

j=2 j=m 



m n 



, fa-(l,^(»,i-l),J,"^) p >^ , 

+ 1^1^ ^2 ■Aim ^i{Pii,j-i),P{i,j))- (3.37) 



i=2 j=m+l ■^1"^ 



- 12 - 



Here we have introduced the shorthand notation 

tr_(a6cd) = (a 6) [6c] (cd) [da] 

and the auxihary functions, 

tr_(m,z,j, l)tr_(m, 1) 



C = 



„2 2 

ti--(l,«, j,m 



(3.38) 

(3.39) 
(3.40) 



Note that b^^^ is symmetric Tinder both i ^ j and 1 <-> m, while v4f^^ is antisymmetric 
under 1 m. The function F^p'' is the finite pieces of the two mass easy box function (or 
the finite pieces of the one mass box function in the hmit where one of the massive legs 
becomes massless). We define the triangle function rj(P, Q) as, 



Ti{P,Q) = Li{P,Q) 



log(P7Q2^ 



p2^o, qVo. 



(3.41) 



(P2 - Q2)i 

If one of the invariants becomes massless then the triangle function becomes the divergent 
function, 



r,(p,Q)^(-i)*- 



e (P2)i 



(3.42) 



The terms associated with a scalar circulating in the loop have the form, 

n m—1 

i=m+l j=2 
m—1 n 

+ EE 



1=2 j=m 



tr-(l,P(ij),i,m)3 



^^4 ^lm-'2(P(j+l,i), P(ij) j 



+ 



tr-(l,P(ij),i,m 



'Im 



1=2 j=m+l 



tr-(l,P(,,,_i),j,m)3 j^i-i) 

+ ^Im ^3[P{i,j-l),P(i,j)) 



tr_(l,P(^j_i),j,m)^^^(^_i) 



/c{^-^^r2(P(,,,._i),P(,,,.)) 



tr_(l,P(,,,_i),j,m) 

-4 Am ^l(-^(jj-l)>AM)) 



(3.43) 



Here we have introduced two further auxiliary functions which are defined as follows, 

'tr_(l,i, j, mf 



-(i^i + i) , 



-jij _ ( tr_(l,i, j, m)^tr_(l, j,i,m) ■ , i 



(3.44) 
(3.45) 
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Figure 3: The MHV diagrams contributing to one-loop i/i-MHV amplitudes. 



3.4 (/i-amplitudes 



The graphs contributing to one-loop (p-MHY ampUtudes are shown in Fig. ^. Diagrams 
of type (b) are the QCD graphs dressed with an additional (p, which may couple at either 
the left or right vertex. The presence of the (p does not alter the spinor structure of the 
amplitudes, so these graphs are exactly those for the pure-QCD amplitudes of the previous 
section, modified to account for the momentum carried by the (p. The ranges of summations 
correspond to those given in eq. ( 3.20 ). 

On the other hand, the diagrams shown in Fig. |^(a) have no counterpart in pure-QCD. 
They all vanish in the limit where the four-momentum of the (p vanishes. The diagram 
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Figure 6: A (/) only diagram in the S2,n channel 



contributing to a cut in the Sj+i^j channel is shown in Fig. ^ 

, . ^ {hi2?{lm)^ ^ ,(0) {i{i + l)){h^2) 

= 4°) + + (3.46) 



with defined in eq. (Q. The dia gram contributing to a cut in the Si+2,i channel 

is shown in Fig. |5|. 

, ^ _ ihmim)' 



+ + l)£2){i2h){h{i + 2)) • • • {iil){hi2) 

= A^^'^ g{i,i + l,i + l,i + 2). (3.47) 

The diag ram contributing to a cut in the S2,n channel is shown in Fig. ^. There are 
contributions from both gluon and fermion loops, and we find, 

(ALAn)^^'' = (§{n, 1, 1, 2) - 4 (l - ^) ^{n, 1, 1, 2) - 2 (l - ^) S{n, 1, 1, 2)) . 

(3.48) 

The diagram contributing to a cut in the Sm+i,m-i channel is shown in Fig. 0. There are 
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Figure 7: A <f> only diagram in the Sm+i,m-i channel 



contributions from both gluon and fermion loops, and we find, 



=A<fi)(^g{m,m + l,m-l,m)-A(l-^^ ^(m, m + 1, m - 1, 



m) 
(3.49) 



Combining all of the diagrams together we find that the cut-constructible pieces of the 
general (j!)-MHV amplitude is given by. 



Cn;i{4'A ,2+,...,m ,...,n+) 

= cr4°) (At^M. ") - - ^) n) - 2 (^1 - ^) 4f,(m, n)) , (3.50) 



where 



n n+i—l 



1=1 j=i+3 



1=1 



(3.51) 



i=l 



We notice that ^^.-|^(m, n) is independent of the position of the two negative hclicity gluons; 
this is exactly as one would expect from an = 4 theory. Nevertheless, the presence of 
the colourless scalar has removed the supersymmetry and as a result we see the appearance 
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of Fg™ functions which are not present in eq. ( 3.36 ). We can write the fermionic pieces as, 



m,ni 



'Im 



--^ml^l(^{j+lj)'^{*.i)) 



m—l n 

^lmF4F''(sij,Si+ij_i;Sj+ij,Sjj_i) 

i=2 j=m+l 
m—l n 

+ ^lmF4F^(Sj,j,Sj+i,i_i;Sj+i,i,Sj,j_i) 
i=2 j=m+l 

2=2 j=m 
m-1 1 

+ E E 

1=2 j=m+l 
m n 

+E E 

1=2 j=m+l 
m—l n 

-E E 

i=l j=m+l 



tr„(l, P(,.,),i,m) ,(^.^1) 

-^im ^UMi,j-i)'Mj,«)-' 

■•^im^l(-P(j+l,j)!-P(j,i)) 





s2 




tr_ 


("1, J) 


,J,1) 








tr_ 




j, m) 





(3.52) 



where the functions 6^^^ and .A^^ are the same auxihary functions as in the pure-glue case 
and are given by eqs. ( 3.39| ) and ( |3.4C1| ) respectively. 
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Finally the scalar pieces are given by, 

m—l n 

4(j>S 



i=2 j=m+l 
n m—l 

i=m+l j=2 



i=2 j=m 



-•^ml-^3(-P(i+l,j).-P(jj)) 



■^ml-^2(^'{j+lj),^'{j,j)) 



+ 



tr_(m,P(jj),i,l)^ 
tr_(m,P(,,,),i,l) ,. 



m—l n 

+ E E 

1=1 j=m+l 



fa-(i,-P(j,i),i,H' 

tr_(l,P(j^^),i, m)^ 



■•^im-^3(-P(j+l,i)'-P(ji)) 



-^lmL2{P{j+l,i), P{j,i)) 



tr_(l, P(j j), J, m) 



Im 



m n 

+E E 

1=2 j=m+l 



+ 



tr_(m,P(^^j),j,l)^^j.(^_i) 



2s4 



tr_(m,P(ij),i,l) j(i-i) 



Im 



^ml ■^2(-P(ij-l),-P(ij)) 



m-1 1 

+ E E 

i=2 j=m+l 



+ 



■^Im L^{P{j,i-l)iP{j,i)) 

tr_(l,P(j^i),i,m)^^^(j_i) 



2s4 



-^Im L2{P{j,i-l),P{j,i)) 



'Im 



(3.53) 



where the auxiliary functions /C^^ and are the same as in the pure-glue case and are 
given by eqs. ( p. 44] ) and (3.45) respectively. 

The similarities and differences between the gluonic MHV and the i?i>-MHV calculation 
are now most obvious. It is clear that both have the same type of auxiliary functions 
multiplying the one- loop basis functions, however the presence of the scalar has introduced 
a second set of summations. One difference is that in the 6-MHV result there are no 
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degenerate triangles. This is a consequence of the absence of C(e ^) terms as predicted by 
the infrared pole structure. 



3.5 Cross Check: The adjacent minus amplitude 

The one-loop (cf), 1~ ,2~ . . . n+) amplitude has been calculated [20] and provides a check of 
our calculation. As mentioned earlier Af^^l is independent of m so we only need explicitly 
check the remaining two contributions, which collapse to, 

A<t>F^r, ^^ _ fa_(l,P(j+i^^),z,2) 
j=3 



and 



n 



+ 2^ — M(-f^[2,i-i),-r(2,i)), (3.54) 

i=4 



tr_(2,P(3,-i),M)% \ 

-^2(P(2,i-l),P(2,j)) I 

, ^ Ar_(l,P(m,n),^,2)3 

tr_(l,P/j, 1 „^,^,2)2 \ 
+ L2(P(m,i),%i)) , (3.55) 

1 



i=3 



'12 



respectively, and which is in agreement with the result of [20] . 



3.6 Cut-completion Terms 



The basis-set of logarithmic functions in which the results are expressed contains unphysical 
singularities, which we remove by adding in rational pieces, the so-called cut completion 
terms. The new basis is given by the transformation, 

Li{s,t) = Li{s,t), 

i.(.,t) = M»,t) + ^(^ + J), 

L3(.,t)=i3(»,t) + ^^(i + ;). (3.56) 

From the breakdown of our amplitude it is clear that only the scalar pieces contribute. 
When considering the overlap terms in the next section it proves most convenient to write 
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the cut-completion terms in the following form, 



CRn{(l),l ,...,m ,...,n+) = r„ 

m n / , ^ \ m— 1 n 



in It / 1 1 \ — / 1 1 \ 

m— 1 n+1 ^ ^ ^ \ m— 1 n / i i \ 

+ E E ^i^'(^a.-i))U:^ + i^ -E E ^i-(^o+M)) jT^ + i^ 



(3.57) 



The factor r„ is given by, 



and 



with 



,a6 _ ("io) (^1) 



(a 6) 

We have also introduced the short-hand notation. 



Ki = - (6 - 6 + 1), (3.60) 

i^^\ = ^^^^^4^-(5-6+l). (3.61) 



Np = 2{l-^]. (3.62) 



4. The Rational Pieces 

In addition to the cut-constructible terms calculated in the previous section, one-loop am- 
plitudes in non-supersymmetric theories also contain rational terms with no discontinuities. 
By definition this means that these terms can only contain simple poles in physical invari- 
ants, which makes these terms amenable to the BCFW recursion relation techniques. So 
far successfTil applications have included amplitudes in QCD [14-16] and the finite and 
adjacent minus 4> amplitudes [50]. 

In an earlier section, we cancelled unphysical poles in Cn by introducing the cut- 
completion terms CRn- If we naively set up the recursion relations we would double count 
on these rational pieces. To avoid this, we define the recursion relation as a function of the 
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physical poles Rn = Rn — CRn- We make a complex shift of the two negative gluons such 
that 

|i) = |1) + 2;|m), |m] = |m] - 2;|1], (4.1) 
ensuring that overall momentum is conserved since 

K(^) =P'i + ^{^ l7l 1] , P^{z) =p'^-^{m m 1] . (4.2) 

The recursion relation on i?„ is defined through the following integral, 

/ = -i, /-(/?„ - CR^). (4.3) 

Provided that z is chosen such that A{z) Q as z goes to infinity, the integral vanishes. 
The residues of the integrand are fixed by multiparticle factorisation so that the rational 
pieces are given by: 

phys poles Zi 

_^ Af{z)RR{z)+RL{z)Af{z) ^ ^^^^^^^^^ CRnjz) ^^^^ 

I ' I 

The final piece of this equation is called the overlap term. It's calculation is relatively 
simple if the poles are all first order. 

4.1 Recursive terms 

The direct recursive terms are obtained by using the following formula 

„D ^ Af{z)RR{z)+RL{z)Af{z) 

= B2 • ^4.5) 

i i 

For our chosen shift ( |4.1| ), the allowed diagrams are shown in Fig. ^, and the summation 
over these given by eq. ( |4Tl| ). In the sum R is defined as the full rational part of the 
amplitude with fewer than n external legs. Due to our choice of shifts the tree amplitudes 

^ W (i+ , i- , - P(- ) , ( j+ , m- , -P+ ) 

are both zero, (here j £ {2,n,{m it 1)}). These three point amplitudes are hence not 
included in either the diagram or the sum. Other terms that vanish are R2{4>y — i") which 
is required to be zero by angular momentum conservation, and R{j~^ , m~ , P^) since the 
corresponding splitting function has no rational pieces. 

Because the tree amplitudes with fewer than two negative helicities vanish, the one 
requires the one-loop contributions with one negative helicity. These are finite one-loop 
amplitudes and are entirely rational. The finite (j) — h • . . + amplitudes were computed 
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Figure 8: The direct recursive terms contributing to -R„(^, 1 ,m , . ■ . , n+) 

for arbitrary numbers of positive helicity gluons in ref. [50]. As a concrete example, the 
three-gluon amplitude is given by, 

Similarly, the pure QCD — h • . . + amplitudes are given to all orders in ref. [14, 60]. In the 
four gluon case, the result is, 

+ + + _ Np (2 4) [2 4]^ 
^4(1 ,2 ,3 ,4 )-^[i2](23) (34) [41] ^^"^^ 

Finally, there the "homogenous" terms in the recursion which depend on the 0-MHV 
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amplitude with one gluon fewer. The first few ^MHV amplitudes are known, 



R2icf>;l-,2-) = -L^(0)(<^^i-^2-), (4.8) 

Rs{<P;l-,2-,3+) = J^aW(0,1-,2-,3+), (4.9) 



i?3(</';l-,2+,3-) = gi2^(°)(<^,l-,2+,3-). (4.10) 



Combining the various diagrams, we find that recursive terms obey the following rela- 
tion, 

n m—1 I 

i=m j=l J+lj' 
jTt— 1 n 



m—1 n -. 

+ E E ^(1"' ■ ■ ■ -^(l+ij)' + 1)^)7^^^°^ (<^' + 1)^' , ■ ■ ■ ,rn-, . . . , z+) 

m—1 n ^ 
j=l i=m *+^'-' 

+R{4>, i-, . . . , (m - 1)+, -P+ (m + 2)+)-J— A(0)(P- , m-, (m + 1)+) 

*m,m+l 

+i?(</>, i-, . . . , (m - 2)+, _,,^), (m + l)+)-^^(°n Vi,-)' ("^ - 1)^'^") 

*m— l,m 

+^(0)(i-, P+ n+)^i?(</<, 2+, , . . . , m-, . . . , (n - 1)+) 

+A(°)(i-,P+2), 2+)-l-i?(0, -P(-2)'3+, , . . . ,m-, . . . ,n+). (4.11) 

The value that z takes is obtained by requiring that the shifted momenta 

^f.-)=^(M)±|H7'^|l]' (4.12) 

is on-shell. In this equation, the sign is positive when the momentum set {pi-,Pj} includes 
pi and is negative when it includes pm- There arc six independent channels, each one 
specified by a particular invariant mass, s^+i^j, Sj^j+i, or by the double invariants, s^.^+i; 
•Sm-i,m) Sn,i and si^2- In each channel, we find that the value of z and the hatted variables 
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are given by, 

Sj+i^i channels 



Zj + l,i 



m) 



|(Pm - -P(j+l,i))^'(j+l,i)l] 



(m|P(j-+i,i)|lJ 



\m 



(m|P(j+i,i)|l] 
Sj,i+i channels z^^j+i : 



\{Pm + -P(j,i+l))P(-,-,i+l)l] 



Sm,m+i channel ^;m,m+i = 

iiv +pm),m + 1] 

1) = ? — ; n , 

' ^ [l(m + l)] ' ' ^ 

^m— l,m chs^nncl •^m— l,m 



[(m + 1) m] 
[(m+l)l] 

= |m + 1]- 



[1 m] 



[1 + 
[(m — 1) m] 

[("^-1)1] 



(m,m+l) 



[1 (m-1)] 



|m — 1] 



[1 m] 



Pi 



(m— l,m) 



5,1^1 channel 
(1 m) 



|1> 



n) 



(n m) 
si^2 channel 
(1 m) 



|1) = |2) 



(2 m) 



\m\ 



Zl,2 



\m\ 



(In) 



(m n) 

+Pm)"-) 

(mn) 
(12) 



(m2) 
|(pi+p^)2) 
(m2) 



P 



[1 (m-1)] 

(mn) 
|2)(mP(i,2)| 



+i)l][m+l| 
[(m + 1) 1] 

|-P(m-l,m)l]["T- - 1| 

[(m - 1) 1] 



p 



(n,l) 



(1,2) 



(m2) 



4.2 The large z behaviour of the completion terms 

In order for the direct recursive contribution to correctly generate the rational terms, the 
shifted amplitude An\z) must vanish as 2: ^ 00. With the shift defined in eq. (4.1) acting 
on two negative helicity gluons this is indeed the case. However, the cut-completion term 
CRn{z) introduced in eq. (|3.57| ) to ensure that the cut constructible part does not have 
any spurious poles, does not vanish as z — > 00. We therefore have to explicitly remove the 
contribution at infinity from the rational part, which now becomes [18, 19], 



Rr, 



InfCP„ 



where 



InfC7P„ = lim CRniz). 



(4.13) 
(4.14) 
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The calculation of Inf Ci?n is straightforward. For the special case of adjacent negative 
helicities, corresponding to m = 2, the cut-completion terms behaves as 1 / z as z — > oo so 
that, 

InfCii:n(</>,l~,2~,...,n+) = 0. (4.15) 

For the general, non-adjacent, case, there is a contribution as z — > oo and we find the 
contribution to be subtracted is, 

Inf Ci?„((;/!), 1~, . . . , m~, . . . , n"'") ^ ^ 



with 



n—l 

Z \'IU ^/ \li llij 11 

m n 



2{m2) {nm)UlZ2{aa + l) . 
1 1 



S,lii"""^''''''^V<m|P,,_,)|l] + <m|P,.)i 1] 



=2 j=m+l 
m— 1 n 



• / 1 1 \ 

S .f+i"''^"'''^'^ V<-|P,„-,)|1] + <m|P,,)|l]J 

m— 1 n—l ^ ^ ^ \ 

^ S£""^'^''''^'^V<m|P,^,,)|l] + <m|P,,)|l] J 

m— 1 n+1 ^ ^ ^ 



i=2 j=m+2 
m— 1 n 



" §.l;i""^'''''^V(-iPo.-i)ii] + <-iP(..)ii] J 

m— 2 n ^ ^ 

^ ?r.Il/"^'^''^^"n<m|Po-+,,)|l] + <m|P(,,)|l]yj' 



=1 j=m+l 



a; 



.o.b 



(m |P(ij) a I w.)^ (am) (6 m) 



2 [la] (a 6) 

and P(^-i) =P(j-i) -pi. 

4.3 Overlap Terms 
The overlap terms are defined as [18, 20], 



(4.16) 



' (^(M)) = ^ .r' w ^ (4-17) 



0„ = ^Res.=.,^^. (4.18) 
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They can be obtained by evaluating the residue of the cut completion term CRn given in 
eq. ( ^.57 ) in each of the physical channels. To expose the coefficients of the poles most 
clearly the cut-completion terms are rewritten as follows, 

m n—1 

CRn = ^ ^ ^ ^ 

j=3 j=m ^'-^ 
n m—l 



m n—1 ^ / 

■_„ •_ ''jj' \ 



m—l 



i=m+l j=2 ^ 

i=3 ^ 

" 1 / •_ 

i=m+l ' 

n-1 . ^ 



i=2 ^^^^ 



' Pml(^{2,n)) - Pml(^(2,n)) 



S2,n 



(4.19) 



with defined in eq. ( |3.59 ). 



The cut-completion terms contain many different simple poles in Sj.j but only those 
invariants which contain either pi or (but not both) have non-trivial overlap terms. We 
observe that the cut completion term contain only simple residues, so for the -P(ij) pole, 
the overlap term is given by, 

(4.20) 



where Zjj- is the value of z that puts -P(j,j) on-shell. The multiplicative factor removes the 
Si^j pole in CRn and replaces it with the correct propagator Sj^. 

The cut-completion terms also contribute to the overlap terms because of singularities 
associated with the multiplicative tree factor in eq. (3.57). The poles in (12) and (nl) 
must be treated carefully, but, since the shift leaves (m| unaltered, there are no overlap 
terms generated by (m (m -|- 1)) or ((m — 1) m). 

Splitting up the cut-completion terms in this way gives the overlap terms the following 
structure, 

m n—1 n m—l m n 

On= EEoUn+ E + E O'n' 

i=3 j=m i=m+l j=2 i=3 i=m+l 

n—1 m—l 

+ E O'n' + E O'n' + O^'" + Of + Of . (4.21) 

j=m j=2 
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We now describe in detail the derivation of each of these terms. 



4.3.1 The overlap term Om,n 

The first overlap terms we consider arc those arising from the Sj j- channel when 3 < i < m 
and m < j < n — 1. Since it is always the case that Pm G -^(ij)) '^^ use the shift 
zi = Sij/ (m |-P(ij) l] • Under this shift the various functions become, 



crNp 



2UIZI (a a + 1) {m |P(,,,) {p, + P(,,,) | 2) <rz | (pi + ) | m> 



, (4.22) 



while, 



<i(^i) 



{ma) {b\{pi+ T\,,j)) P^,^J^J\|-n) 
{ah) <m|P(,,,-)|l] 
{ma)''- {h\{pi + Pj^.^j)) Pt^,,^j)\mf 
(a6)2(m|P(,,,.)|l]' 



6+1) . 



(4.23) 



The prefactor multiplying the A and K functions is simplified since P(i,j) in the numerator 
is never shifted (as it is always adjacent to a (m|). 



(m |P(ij) a\ i)" _ (m \P{i,j)\a] /(a |(pi + P(i,j)) P(i,j)| m)" 



(a| ^ |a]"-i l%,)l 1] V {a |P(ij)| 1 



(4.24) 



Oli^n is thus given by, 

Tn{zi) 



1 ( {j\iPl+P(^,J))Pii,J)\mf j^i^i) 

s^A{m\P,M[ 3 0- IP,,) 1 1]^ ^''^ 

+ 2<j|P(,,)|l] j + - .7 + 1, ^(.J) - PiiJ)) 

{m \P(i,j)\i] ( {i \{pi + P(,,)) m)^ 

(^|(pi + ^(»,i))^(»,i)|m)' . . A . . 

— ^mU^l) I + ^ « - P{i,j)) 



2<^I%.)|1] 



(4.25) 



4.3.2 The overlap terms O^", On'^ and O^'" 



The contributions in the Si^n, S2j and S2,n channels are evaluated under the same shift as 
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Om,m such that, 



OXn= rn(^l) 



O: 



2,n 



1 r(m|P(j_„) 



.Sj,„l(m|P(j,„)|l] 



(^1 \P(i,i) Pii,n)\m)^ {mi) (n|P(j,i) P(j_„)|m) 



3(i|%„)|l] 



(in) (m|P(j,„)| 1] 



2(i|P(,,„)|l] (^n)'<m|P(,,„)|l]' 



P(i,j)) 



(m|P(i,„)|l] V 3(n|P(,,„)|l] 



A 

2 ^ml 



\2 

/ r^n(i-l) 



l^l) + o/ IP Ml -'^"^l 



(^i) 



1 

S2j 



{m \P(2,j)\j] ( <i l^(i,i) P{2,j)\mf {mj) (2 |P(i,j.) P(2j)| m> 



[m\P( 



(2,j) 



3<il^(2,,)|l] 



<J I^'CU) P{2,j)\ m) {mjf (2 |P(ij) P(2j)| m)' 



2<il^(2,,)|l] 



(i2)2(m|P(2,,-)|l]' 



(j2) <m|P(2,,-)|l] 
+ {j^j + hP{i,j)^P(i,j)) 



{m |P( 



(2j) 



<"^I^(2J)|1] 



ii,j)P{2,p\m 



3<2|P( 



(2,i) 



2(2|^'(2,,)|1J 



rn(^;i) 



1 r(m|P(2,„)|n] / (n|P(i^„)P(2,„)|m)^ (mn)(2|P(i^„)P(2,„)|m) 



(2n) <m|P(2,„)|l] 



+- 



Ls2,nl<m|P(2,„)|l] V 3<n|P(2,„)|l] 

{n \P{i,n) P{2,n) I m)^ (m n)^ (2 |P(i^^) P(2,n) I m) 
2<n|P(2,„)|l] (2n)2<m|P(2,„)|l]' 

(m |P(2,,ol 2] / (2 \P(i,,) ^(2,,ol '">' 2) (" 1-^(1.0 -^(2,.)l 



<m|P(2,„)|l] V 3<2|P(2,„)|1]' (2n)<m|P(2,„)|l] 

(2 |P(i,n) P(2,n)l (m2)^ (n |P(i,,) P(2,n)| m)^ 

2<2|P(2,„)|1] (2n)2<m|P(2,„)|l]' 



(4.26) 



4.3.3 The overlap terms O^'^, On'^ and Oj;^ 

A similar set of overlap terms are generated in the Sij, Sij and Si i channels when pi € 
P(ij)- We therefore use the shift Z2 = — •Sij/ (m|P(jj)| l] . Once again the tree factor, F 
and the functions A and K must be evaluated under this shift; 



^n{z2) 



crNp 



, m 



with. 



2n^-^ (aa + 1) <m |P(,,,-) (pi - P(,,,))| 2) {n\{pi - P(,,,)) P(,j)| m) ' 



Atmiz2) 



(6m)(m|P(.^j)(pi-P(^j))|a) 
(a6) (m|P(„.)|l] 

(6 m) (m| P(jj) (pi - P(ij))| g) 
(a6)2(m|P(,,,-)|l]' 



(6^6+1) , 



(6^6+1) 



(4.27) 



(4.28) 
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Finally the prefactor multiplying the A and K functions is given by, 



(ljP(,,,)a|mr 



(-1) 



(a|P(,,,)|a]"-i ' ' {m\Pii,j)\l]\ (a |P(,,,)| l]""' 
The overlap contributions are given by, 



(4.29) 



, 2 ^Im 



<"^l%i)ll] 



"^l%i)|lJV 3(j|P(,,,-)|l] 

KZ~"i^2)]+{j^j+l,Pii,j)^Pii,j)) 
2 ^lm(^2) 



-Pi)^ {jmf 
2 0-|%,)|l] 



{P(i,j) -Pi)^ (^"i)^ 



^l4(^2) +(^^^-l,P(^J)^P(^J)) 



On'' = r„(z2) 



(m |P(i,j)| n] / (P(2,j))^ (nm)' 



sijl3(m|P(2j)|l] V (n|P(2j)|l 



^l4(^2) 



(P(2.j))^ (nm)- 
2 1^(2,) 1 1] 



a,3 



+ - 



{P(2j)f {jmf /{nm) (m - P„j) P(i,,)| j> 



3 0-|^(2,,)|l] 



(jn) <m|P(i,,-)|l] 



+ (Im) 



, (-P(2,j))'^ 0'"T-)^ {nmf {m\{pi - Pnj)P(ij)\jf 2 

+ r,, ,9 



2<J|^(2,,)|1] V (in)^(m|P(i,,)|l]' 



m / {P{i,n)f {irnf / (2 m) {m |P(i^„) P(j,i)| i> 



,il <m|P(,i)|l] V3(^|P(,„)|1]' V (2i)<m|P(,,i)|l] 
iP{i,n))'^ {imf f {2mf (m|P(i^„) P(i,i)\if 



+ (Im) 



2(^|i='(^,n)|lJ 



(2i)^(m|P(,,i)| 



+ (im)^)) + (i ^ z - 1,P,,„ ^ P,,„) 



(m|P(,,i)|2] / (P(,,„))6(2m)= 



2(i-l) 



3(2|P(.,n)|l 



_ A- 
2 ^Im 



(Z4) + 



(P(,,„))" (2 m)- 
2<2|P(.,n)|l] 



if 



2(i-l) 
Im 



(^4) 



(4.30) 



At first glance there appear to be poles of order greater than one in S12 and however 
the presence of the factor {P(ij) — Pi)^" ensures that when i = 1, j = 2 or i, j = 1 there 
are no issues with higher poles. As a result, poles in these channels are only generated by 
the multiplicative tree factor. 



4.3.4 The overlap terms O^" and O^^ 

The final two overlap terms have a more subtle origin than the previous contributions 
and since they come from the tree factor the form of the cut-completion terms, can be 
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compacted. Wc will however, need to have forms for Sij when they acquire a z dependence 
(with z now in either the S12 or s„i channel). We consider first the S12 channel (2:3 = 
— (1 2) / (m2)), we use the following form for Sij, when pm € P{ij)' 



Si J (^3) = 



{m\P(i,j) (P(.,,)+pi)|2) 
(m2) 



We will also require the tree factor T and the functions A and K: 

crNp (2 m) 



2n^-i {aa + 1) (Im) (n2) (12) 



( {h2){lm){ma) 
V {2m) {ah) ^ ' 
(62)^l^)^^a)^ 

(2m)2(a6)2 ^ ^ 



(4.31) 

(4.32) 
(4.33) 
(4.34) 



The overlap terms associated with this channel are defined by using the (1 m) symmetry 
of the cut-completion terms: 



(^12 _ (^12 I (Qii 



12 

n 



(4.35) 



With 0^2„ defined by, 

m n 
i=3 j=m 



{m\P^i,j)j\2f{lm)' 



i(i-i) 



3 (m \j{P^i,j-i) + Pi) + P(i,j-i)j\ 2) 



-A- 



(^3) 



+ 



H%,)i|2)'(lm)^ 



2 (m |i(P(i,,-i) +Pi) + Pii,j-i)j\ 2) 
1 



.m 



\P(:,j) (P(„-)+Pi)|2) 



+ 



;^r'^(-3) 
1 



l^j-i) (^(M-i)+Pi)|2> 



m n— 1 

+ 1:1: 

1=3 j=m 



(m|P(,,,.)i| 2)^(1 



m) 



3 (m K(P(i+ij) +Pi) + P{i+i,j)A 2> 



r^^i(^3) 



<m|P(,,,)z|2>'(lm)^ 



2 (m K(P(i+i,j) +Pi) + P{i+i,j)i\ 2) 
1 



i^i(^3) 



+E 

1=3 



{m \P(ij) {P(ij) +pi)\ 2) (m (-P(i+ij) +Pi)| 2) 

(m|P(i,„)i|2>'(lm)=^ (mi)(n2)(lm) 



3 (m KP(,+i,i) + 2)' n) (2 m) 

(m i\ 2f (1 m)' (m i)' (n 2)' (1 m)' 

2 (m + P(i+i,n)«l 2) (i (2 m)' 

1 1 



+ 



(m I 2) (m |P(i+i,n) -P(i+i,i) I 2) 



(4.36) 
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For the second set of sums we will need to know Sij with pi G P{i,j) , 



(m2) 



(4.37) 



We also require, 



{bm) (Im) (2a) 
(2 m) (a 6) 
{bmf {Imf {2af 
{2mf{abf 



-{b^b+1) 



{b^b + 1) 



(4.38) 
(4.39) 



so that, 



^12 



r„(^3) 



n m—1 



i=m+l j=2 



{2\P^,j^j\m)'{lmf 



3 (m |i(P(i - Pi) + P^i,j-i)j\ 2) 



-A- 



(^3) 



(2|P(,,,-)i|m>'(lm)^ 



2 (m I j(P(M-i) - Pi) + P(*j-i)j| 2) 



+ 



K 



1 



<m |P(,,,.) (P(,,,.) - pi)| 2) (m |P(,,,._i) (P(,,,._i) - pi)| 2) 
I y^' f (2|P(i,,.)i|m)^lm)=^ / (2 j) (nm) (1 m) ^ \ 

,t^l 3<m|iP(2,,_i)+P(i,,_i)i|2>'V {jn){2m) ') 



3<m|iP(2j_i)+P(i,,_i)i|2) 

(2| P(i^^.) il m)' (1 m)^ / (2 (n m)" (1 m)- 



2 (m I jP(2j-i) + P(i,,-i)j| 2) V (j n)2 (2 



m) 



(1 m)^ 



+ 



(m |P(i,,) P(3,,)| 2) (m |P(2,,-_i) P(3,,--i)| 2) 



n m—1 



+ E E 



v2|P(^,/) 'I "'>'^ (-'-"') 



m+l j=2 

I 3 (m K(P(i+ij) - Pi) + P{i+i,j)i\ 2) 

{2\P{i,j)i\mf {^mf 



2 {m \i{P(i+ij) - pi) + P{i+i,j)i\ 2) 
1 



+ 



{m |P(ij) (P(ij) - pi)\ 2) {m \P(i+ij) {P{i+i,j) - Pi)| 2) y . 



(4.40) 



In writing the above, we have used ^^^(23) = 0. 

The final overlap term is O"^ and is calculated noting the (n ^ 2) symmetry in the 
shift. Once again we define the usual functions under the shift (z^) = — (Ira) / (mn) , 
When pm e P^ij), 



[mn) 



(4.41) 
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together with, 

r(z,)- (4 42) 

.r,h , N f (b n) (1 m) (m a) , x\ , ^ 

^^(.4) =( ""'f<\^f (4.44) 

\ (nm) {ao) J 
The overlap terms in this channel are again split into two terms 

Of = Ot,n + Oil, (4.45) 

with 

QUI _ V Vr {zA[ (m|P(,,,-)i|n)=^(lm)^ A^^'-^\~^ 

2(m|j(P(.,,-i)+Pi) + %,-i)j|n> 



1 1 

+ 



1%^-) (P(i,^-) 2) (m +Pi)| n> 

, y [ H^(2,,)i|n)'(lm)^ {mi)(2n)(lm) 
.t;;, l3<m|iP(i,,_i) +P(2,,_i)i|n>' (i2) (nm) 

(m| P(2,,) il n)' (1 m)^ (m jf (2 n)^ (1 m)^ 
2 (m |jP(ij-i) + P(2j-i)i| n) (j 2)2 (n m) 

1 1 

+ 



2 



(m |P(2j) P(ij) I 2) (m |P(2,j-i) -P(i,j-i) I n> 

m n— 1 /In -I /i \3 

(m|P(jj)z|n) (Im) 



1=2 j=m 



3 (m K(P(i+i J) + Vi) + P{i+i,j)A ny 



2 (m K(P(j+i J) + Pi) + P(i+i 



j) 



z n 



-K!l(^3) 



+ 



(m |P(i J) (P(i J) + ^?i) I n) {m |P(i+i j) (P(i+i,j) + Pi) I n) 



For the second set of terms we need to evaluate Sjj when pi G P(ij), so that. 



and. 



(4.46) 



{m n) 



<fa)= ( '"7"'";""°' -'"-"+'>)• <^-'«" 

g&W=( """>'",rf <"°'' -(6^t+l)). (4.49) 
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We find that, 



Ol,n= r„(z4) 



n m—1 

E E 

i=m+l j=2 



{n\ P{i,j)j\ rnf {ImY 



3 (m I - Pi) + P{i,j-i)3\ n) 



— A- 

2 Im 



(^4) 



{n\ P{i,j)j\ m)^ 



2 (m I - Pi) + P{i,j-i)j\ n) 

1 



K- 



+ 



1 



{m \P(ij) (P^ij) - pi)\ n) (m \P(i,j-i) {P{i,j-i) - Pi)\ n) 

n m—1 ^ /In -I \3 /1 \3 



+ E E 



i=m+i j=2 ^ 3 (m - Pi) + Pii+i,j)i\ n) 

(n \P{i,j) i\ rnf' (1 m)'^ 



2 (m - Pi) + P{i+i,j)i\ n) 

1 



^rm(^4) 



+ 



+ E 



|-P(i,i) ^1 m-)^ (1 rn)^ 



i=m+l 



^^^^ ^ (ni) (Im) (2m) 
3 (m + n)' V 2) (n m) 



(n|P(,,i)i|m> (lm)2 / (n i)^ (1 m)^ (2 m)^ 

' (Im) — 



2(m|iP(j+i^„) 



71 / 



(i 2)^ (nm)^ 



1 



+ 



1 



("i l^(i,l) ("1 l^(i+l,l) P{i+l,n)\'n) J _ 



(4.50) 



5. The four point amplitude 



The calculation of all Higgs plus four-gluon amplitudes at NLO in the heavy-top effective 
theory has been performed numerically in [61]. Here we provide an analytic form for the 
A^I\h,1-,2+,3-,A~^) to illustrate the use of our results for the (p-MHV amplitude for 
general n. 

The cut-constructible part of the cp-MHV four point amplitude is given by setting 



n = 4 and m = 3 in eq. (3.50), using the gluonic, fermionic and scalar contributions given 
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in eqs. (|3.51 ), ( 3.52 ) and ( 3.53| ) respectively, 



f 1 

CiifpA 5 2"*", 3 ,4:'^) = < - - ^ F|™^(si,i+3, Sj+i,i+2; Si+i,i+3, Sj,j+2) 

^ i=i 

1 ^ ^ 



-4 1 



Nf 



ltr_ (3241) tr_ (3421) 

r 4F IS234 , S23 , ■S34 J 



2 s2 s2 



tr_ (3241) tr_ (3421) 



S24S?3 



Liis23, S234) + (2 ^ 4) + (1 ^ 3) + (1 ^ 3, 2 ^ 4) 



-2 1 



Nf 



ltr_ (3241)2 tr_ (3421)2 



4 ^4 



Sn/i S 



F4F (S234, ■S23, ■S34) 



24 •'13 



tr_(3241)tr„(3421) /tr_(3241)2 

:4 I — -t^3(S23, S234) 



'13 



+ 



tr_ (3241) tr_ (3421) 



2s2 

^*24 



^2(5235 ■5234) 



3S24 

tr_ (3421) tr„ (3241) 



-^l('S23) -5234, 



'24 



+(2 ^ 4) + (1 ^ 3) + (1 ^ 3, 2 ^ 4) 



The cut completion terms are given by eq. (|3.57] ) , 



(5.1) 



Ci?4(</',l-,2+,3-,4+) 



Np 



1 



327r2 (12) (2 3) (34) (41 
(3|24|1)^ (34) (21) 



(3|24|1)^ (34)^(21)2 



3(S234 -S23)2 (4 2) 2(S234 -S23) (4 2)^ 

+(2 ^ 4) + (1 ^ 3) + (1 ^ 3, 2 ^ 4). 



1 1 

+ 



S23 ■S234 



(5.2) 



The remaining rational contributions are obtained by shifting the two negative helicity 
gluons, 



|i) = |l)+z|3), |3] = |3]-z|l]. 



(5.3) 



As discussed in subsection this shift generates a non-vanishing contribution as z — > oo 



in the cut completion term CR4. To compute this contribution, we use eq. ( 4.16 ) with 
m = 3 and n = 4 to find. 



InfCi?4(</',l",2+,3",4^ 



Np (2 3) (3 4) [2 4] 2 
'327r2 (2 4)=^ [12] [41] 



(5.4) 



The direct rational contribution is generated by the recursion relation ( 4.11| ), again 



-34- 



with m = 3 and n = 4 and is given by, 

/?4(</',l",2+,3-,4+) = A(o)(0,i-,p-34)— i?(-P2+34,2+,3-,4+) 

■5234 



+i? (</., i- , 2+ , -P3+4 , ) — ^ W (P3-4 , 3- , 4+ ) 

i~ , 4+ , -P2+3 ) — ^(0) (P2-3 , 2+ , 3- ) 

+yl(°)(i-,P+,4+)— i2(,^,-P4l,2+,3-) 
S41 

+A(°)(i-,P+>2+)— i?(,/),-Pr2,3-,4+), (5.5) 

•512 

where we recycle the known lower point amplitudes. For the four-point amplitude, we 
require the rational parts of the <f> with one minus and two positive helicity gluons ([4.6|), 
the two and three-point 0-MHV amplitudes given in eqs. (|4.8|) , ([4.9|) and ( 4.1[l| ), as well as 
the pure four-gluon QCD amplitude with a single negative helicity of eq. 

We find that 



Similarly, 



^234 _ Np m% (24)[24](3|P234|1]' 

" 967r2s234 (4|P234|l]'(2|P234|l]'' 



Rf = Rf (2^4). (5.7) 
In the other channels, 

i?f = -^AW(cA,l-,3-,2+,4+) 

Rf = Rf (4^2), (5.8) 

and finally, 

o4i2_ Np [2Af (3IP412II]' 

" 967r2 S412 (2 4)[12]2[41]2- ^"'-'^ 

The overlap terms are given by, 

04(</., 1-, 2+, 3-, 4+) = Of 4 + Of + Of + Of + Of + Of\ (5.10) 
The first term is generated by eq. ( [4.26 ) with n = 4 and has the following form 

Np /l(3|P234 Pl234|2)2[4 2] 



^ ~ 327r2s234 V3 (2 4)(2|P234|1]' 

l (32)(3|P234 Pl234|2)(3|P234 Pl234|4)[42] ^ ^ ^^ | 

2 (2 4)2(2|P234|1](3|P234|1] 
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The overlap pieces in the 23 and 34 channels are given by eq. (|^) and eq. (|]2|) (with 
« = j = 3), 

q23^ Np ( (32)^(4|Pi23|2]^[24] 
^ 327r2s23V 3(4|Pi23|l]^(4 2) 

(32) (34) [24](2|Pi23|2](4|Pi23|2] \ 
2[12](42)2(4|Pi23|l] 

034= q23 (4^2). (5.13) 

0^1 and both vanish, while eq. (4.30) with i = A, j = 2 leads to, 



or = -^('-^''^ ^''^ + (2^4)). (5.14) 



327r2s4i2V2 (2 4)(3|P4i2|l][41] 3 (24) [41] 

Combining contributions, the full four-point amplitude is given by, 

4'^ (0, 1- , 2+ , 3- , 4+ ) = O4 (0, 1- , 2+ , 3- , 4+ ) + CR^i^, 1" , 2+ , 3" , 4+) 

+ i?4(0,l^,2+,3-,4+), (5.15) 

with 

%,l~,2+,3-,4+) = 04(</.,l-,2+,3-,4+) + i?4(</',l~,2+,3-,4+) 

- InfCi?4(</',l~,2+,3~,4+). (5.16) 

After some algebra, the combination of overlapping and recursive terms can be written in 
the following form, free of spurious singularities ^, 

ii4(</',l",2+,3-,4+) = -^A(0)(Al~,2+,3-,4+) 

Np [2 4]"^ / _ 523^34 o£23£34 _ 512^41 o£l2£41 \ 

1927r2[12][2 3][34][41]V S24S412 sj^ S24S234 sj^ /' 

(5.17) 

where A^*^) (A, 1^, 2^", 3^, 4^) is the difference of and i;^^ amplitudes. Finally the full Higgs 
amplitude is given by the sum of (p and cj)^ amplitudes 

4')(i/,l",2+,3-,4+) = 4^)(<A,l-,2+,3-,4+) + 4^)(<A^l-,2+,3-,4+), (5.18) 

with, 

4')(</)^l-,2+,3-,4+) = 4^)(</),2-,3+,4-,l+)(,,)^[,,-]. (5.19) 



We note that the rational terms not proportional to Np in eq. (|5.17 ) cancel when forming 
the Higgs amplitude, just as for the 4^'*(f/', 1^, 2^, 3"^, 4+) amplitude of ref. [20]. 



^Which we have checked with the aid of the package S@M [62] 
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The one-loop amplitudes in this paper are computed in the four-dimensional hehcity 
scheme and are not renormalised. To perform an MS renormahsation, one should subtract 
an MS counterterm from 

A^^^A^)-crf-^Af. (5.20) 
The Wilson coefficient ( p.2| ) produces an additional finite contribution, 

4^)-4^) + y 4°). (5.21) 



6. Cross Checks and Limits 

6.1 Infrared poles 

The infrared pole structure of a one-loop (?!>-amplitude has the following form, 

4i) = _^^(0)J-r_i^^ +0{ey (6.1) 

Since only A'^^{m,n) contributes at 0{e~'^) the IR pole structure of the general (/)-MHV 
amplitude is identical to that of the adjacent minus case (apart from the trivial change in 
the tree amplitude). This combination was shown to have the correct IR behaviour in [20]. 

6.2 Collinear limits 

The general behaviour of a one-loop amplitude when gluons i and j become collinear, such 
that Pi — > zK and Pi+i ^ (1 — z)K, is well known, 

4^)(...,.^',. + i^-s...)^^^ 

[ (. . . , i - l^'-i ,K^,i + 2^'+2 , . . . )Split(°) (-i^-^; , i + l^»+i ) 

h=± 

+A^^l^{. l^'-i ,K'',i + 2^>+2 , . . . )Split(^) i-K-''- + l^^+i )] . (6.2) 

The universal splitting functions are given by [1, 2, 63], 

Split(0)(-i^+;l-,2+) = , (6.3) 

7^(1^(12)' 

Split(°)(-i^+;l+,2-) = . I] ~ , (6.4) 

^z(l - z) (12) 

SplitW(-i^-;l+,2+) = (6.5) 

V^z(l-2)(12) 

Split(°)(-i^-;l-,2-) = 0. (6.6) 
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The one-loop splitting function can be written in terms of cut-constructible and rational 
components, 

Sp\it'^^\-K-^,1^\2^^) = Split(^)'^(-K-'*,l^i,2^2) + Split(i),i?(_^-ft^lAi^2^2) j^gj) 
where 

Split«'^(-i^±,l-,2+) = Split(°)(-i^±,l-,2+)^ X 

(i - , (i, 1 ) - (i. 1 - ^))) . (6.8) 

SplitW'^(-i^+,l-,2-) = Split(°)(-i^+,l-,2-)2L X 

Split(^)'^(-i^-,l-,2-) = 0, (6.10) 
Split(i)'^(-i^±, 1", 2+) = 0, (6.11) 

SpU.«.'=,-i.M-,2-) = ^^^, (6.12, 

SpU.<-'^''(-i.-,l-,2-)^^ V^(12) (,^3, 

Explicitly, the cut-constructible parts should satisfy, 

C„(...,^^^^ + l^-^...)^^^^ 

/l=± 

Cn_i (. l^^-i , K'', i + 2^'+2 , . . . ) Split(°) {-K-''; + l^'+i ) 

+aJ2i(. . . , i - l^*-i , i^'', i + 2^>+2 , . . . ) Split(^)'^(-ii:-''; i^' , i + l^'+i), (6.14) 

while the rational pieces obey, 

h=± 

Rn-i{. 1^'-' ,K^,i + 2^'+2 , . . . ) Split(°) {-K'^; + 1^'+^ ) 

+aJ2i(. . . ,z - \^i-\K^,i + 2^'+\ ...) SplitW.«(-i^-'^;i^i,i + i^i+i). (6.15) 

6.3 Collinear factorisation of the cut-constructible contributions 

In Ref. [20], it was demonstrated that the helicity independent cut-constructible gluonic 
contribution obeys, 

h=± 

C^l? (• . . , i - 1^'-^ ,K'^,i + 2^^+' , . . . ) Split(°) i-K-^] + 1^'+^ ) 

(. . . , i - l^'-i ,K'^,i + 2^'+2 , . . . ) Split(^)'^(-ii:-''; + l^'+i ). (6.16) 
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Therefore to check the coUinear behaviour of the general i;^-MHV amphtude, we simply 
need to check that the fermionic and scalar contributions satisfy the following relation, 

cS'^H- . . , i - 1^'-^ ,K'',i + 2^'+^ ,...) Split(°) {-K~^; + 1^'+' ) . (6.17) 

h=± 

In other words, the F and S contributions should factorise onto the tree-level splitting 
amplitude for the helicity of the gluons considered. According to the definition of C„ in 
eq. ( ^ ), there is an overall factor An \ which in the collinear limit produces the correct 
tree-level splitting function. It therefore remains to show that, 

- <-t! (6.18) 

in the collinear limit with Aff^{m, n) and A'^^{m, n) given in eqs. (|]5|) and ( ^1% respec- 
tively. 



6.3.1 Collinear behaviour of mixed helicity gluons 

We first consider the limit where two adjacent gluons become collinear, one of which has 
negative helicity. For definiteness, we take the limit {m — 1) || m. 

The coefficient of the box function l)'^^ enters both A'^^ and A'^'^ . In this limit, 

Aj rn-l\\m tr, (i^, i, j, 1) tr. (JC, j, i, 1) _ ij 

"ml ' 2~~2 — "KV (^D.iyj 

For the special cases, i = m — 1 and j = m — 1, we have, 

Cr'" = eT"' = (6.20) 
so that the box contribution correctly factorises onto the lower point amplitude. 

The remaining terms in the sub-amplitudes are proportional to one of the auxiliary 
functions .F^^ with F = A^K, and I and which are defined in eqs. ( p. 401 ), ( p. 44 ) and ( 3.45| ). 
We shall see that these too have the correct factorisation properties. Let us first consider 
the ranges 2 < i < m — 1 and m < j < n. When i < ni — 2, the momentum P[ij) always 
contains both m — 1 and m, while P(j^i) never includes either m — 1 or m, and we find 
relations such as, 

tr_(m,P(^j),i, 1) m-i\\m tT_{K,P^ij),i,l) ^-^ 
tr„(l,P(^-^),i,m) m-i||m tr^V^O^^ 

We note that for the special case i = m — 1, 

.m-i,j _ tr_(m,j,m - 1,1) _ tr_(m, j,m, 1) m-i||m 

,Tn — 1 -i m— lllm 



m—l\\m 
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Similar relations hold for the terms involving /C and X. Therefore, all terms in the n-gluon 
version of and therefore either collapse onto similar terms, or vanish in such 

a way that the reduced summation precisely matches onto the corresponding vl^^^.^ and 

a4>s 

"^n-l;l- 

6.3.2 Two positive collinear limit 

Next we consider the limit when two positive helicity gluons become collinear. We focus 
on the specific example where ^ — 1 || £ with 3<£<m — 1. As in the previous subsection, 
let first consider the ranges 2 < ? < m — 1 and m < j < n. We note that, 

bt b^- (6.23) 

The collinear factorisation of box functions has been well studied [1, 2, 63] and in this case, 
the relation, 

{bim^) ^iTist-ij, se,j-i; se,j, se-ij-i) + (b{^^^ Fl^^{sej, s^+ij-i; s^+ij, sej-i) 

* [blm) ^ 4F \SK,j, Se+l,j-l, SK,j, S^+lJ-l j 

(6.24) 

ensures the box terms correctly factorisc onto the lower point amplitude. 

The next set of functions we consider are the triangle functions which have j as the 
second index, these functions possess the general form: 

e 

J2 tr-(m, P(,,,), j, 1)"J-4^')l„(P(„_i), P(„)). (6.25) 
i=e-i 

There is no contribution when i = I, because J^i^\ = J^{^ = 0, while the remaining 
i = i — 1 contribution collapses onto the correct term, 

tr_(m, P^K,j)J - 1, ir^i^i~'^LniP^K,j-i), P(K,j))- (6.26) 

Similarly, when we consider 

i 

J2 tr-(m, P(,- irTilLn{Pij+i,i),P^j,i)), (6.27) 

i=e-i 

there is no contribution when i = £ — 1, while for i = £, we recover the correct contribution. 

The remaining types of triangle function are of the form 

e 

J2 te-(m, P(,,,), z, l)".^^iL„(P(,+i,,.), P(,,,.)). (6.28) 

i=£-l 
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Since J'^i = J'^i have contributions from both terms, however, it is straightforward 

to show that, 



tr_(m, ^ - 1, l)"'LniP(e,j),P{e-i,j)) + tr_(m, P(^+ij),£, l)''Ln{P(e+i,j), P{e,j)) 



tr_(m, P(^(.+ij),K, l)"-Ln{P[e+ij),P(K,j))- 

(6.29) 



Similar considerations apply to 



tr_(l,P(,-,),i,m)'^.F^^ (6.30) 



thus ensuring the correct collinear factorisation. 



6.4 The cancellation of unphysical singularities 



The cut constructible terms eq. ( 3.52| ) - ( 3.53| ) contain poles in (ij). For the most part. 



and j are non-adjacent gluons and as such there should be no singularity as these become 
collinear. In the following section we prove that this is indeed the case. To be explicit, we 
consider the collinear limit i \\ j with. 



i — > zK, 

(1 - z)K. (6.31) 



Let us consider the cut-constructible pieces associated with the fermionic loop contri- 
bution, Aff^{m,n) given in eq. (im). There are ten terms containing an explicit pole in 
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Sij which are given by, 



(m,P(i+ij),i,l) tr_ (m,i,j,l) 



tr_ 



tr_ 



tr_ 



tr_ 



tr_ 



tr_ 



'Im 



m, P{i+i,j-i) , i, 1) tr_ (m, 



l,P(j-i_i),i,m) tr_ {l,i,j,m) 



Li(.P{i+ij-i),P(i,j-i)) 



Li{Pu,t-i),Pu,t)) 



Im 



tr_(l, 



Li{P(j+i,i-i), P(j+i,i)) 



m 



,P{i,j-i)JA) tr_ (m,j,i,l) 



'Im 



'Ti, -P(j+i j-i) , J, 1) tr_ (m, j, i, 1) 



hP{j+i,i),j,m) tr_(l,j,i,m) 



'Im 



1, P{j+i,i-i) > J> "i) tr_ (1, j, i, m) 



'Im 



^l(^'(j+M-l),-PO-,j-l))- 



(6.32) 



Using P(i+ij) = P(i+ij_i) P(j,j-i) = P(j+i,i-i) +Pj, P{ij-i) = P{i+i,j-i) +Pi and 

-P(j+i^i) = -P(j+i^i_i) as well as tr_ (l,j,i,m) = — tr_ {l,i,j,m) + O(sij) etc, we can 
rewrite these terms as 

X - ^l(^(i+lj-l)>%j-l))) 

X (i^l(P(iJ-l),P(i,j)) -i-l(P(i+l,j-l),P(i+lJ))) 

X (Li(P(j+i,j),P(j,i)) --Li(P(j+i^j_i),P(j^j_i))) • 

(6.33) 

Finally, in the i \\ j collinear limit, 

tr_(m,P(j+ij_i),i,l) (Li(P(j_^ij), P(jj)) - ii(P(.j+ij-i), P(i 
^ tr_(m,P(j+ij_i),i,l) (Li(P(jj_i),P(jj)) -Li(P(j+ij_i),P(j+ij))) (6.34) 



tr_ 


(m,P(i+ij 


-i),i,l) tr_ 


("^,^, j,l) 










tr_ 
+ — 


(m,P(j+ij 


-1), j,l) tr_ 


("^,«,i, 1) 








Sij 


tr_ 


(l;-P(i + l,j- 


i),i,m) tr_ 






''Im 




Sij 


tr_ 
+ — 




i),i, m) tr_ 










Sij 
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and noting that the combination, 

(6.35) 

we see that all singularities cancel. The same arguments apply to the cut-constructible 
pieces associated with the scalar pieces. 

6.5 CoUinear factorisation of the rational pieces 

This section is devoted to the collinear factorisation of the rational pieces of the four point 
amplitude. Since there is a (1 3) and (2 •s-^ 4) symmetry there are two independent 
limits 1 II 2 and 2 || 3. We first consider the collinear limit 2 || 3. It is straightforward to 
see that the amplitude correctly factorises onto: 

i?4(<^,l",2+,3-,4+) + Ci24(</',l-,2+,3-,4+)^i^i?3(0,l-,/C+,4+)Split(o)(-K-,2+,3-) 

+i23(0,l-,i^-,4+)Split(o)(-iC+,2+,3-). 

(6.36) 

In a similar fashion the remaining non-trivial collinear limit takes the form, 

^4('^,l-,2+,3-,4+) + Ci?4(<^,l-,2+,3-,4+)'i'i23(./',i^+,3-,4+)Split(0)(-K-,l-,2+) 

+i?3(.^,if-,3-,4+)Split(°)(-i^+,l-,2+). 

(6.37) 

6.6 Soft limit of A^^^((^, 1", 2+, 3", 4+) 

The final test is to take the limit as the (f> momentum becomes soft. Our naive expectation 
is that in this limit, the 4> field is essentially constant so that 

Ccf>tiGsD ^^uG^sD ^ trGsD ^,uG^sD ■ (6-38) 

In other words, the amplitude should collapse onto the gluon-only amplitude. Follow- 
ing [50], we expect that, 

4^) (</., n-g- , n+g+) n_ {n^g' , n+g+), (6.39) 

while 

4)(</)t, "'-^^ n+ A^^\n^g- ,n+g+). (6.40) 

Wc first consider the cut constructiblc contributions. These factorise onto the four 
gluon amplitude in rather trivial manner since in our construction we separated gluon-only 
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like diagrams and those which require a non-vanishing (/)-momentum. In the soft hmit, the 
one and two mass easy box and triangle functions have smooth limits so that, 



2 

-ml 



0, (6.41) 



"'^"O. (6.42) 

(6.43) 

Furthermore, in the soft limit the L/. functions become the massless Tj functions defined 



in eq. (3.42), 



r ( N Bub g234 - Bub g23 p^^o -1^ ( /^^ , ... 

Lk{s234, S23) = 7^ — • (6.44) 

(5234 -523)" sl^e{l - 2e) \- S22, ' 



Altogether, we find that 



C74((A,l-,2+,3-,4+)''C°2C4(l-,2+,3-,4+), (6.45) 



where 6*4(1 ,2^,3 ,4"^) is given by eq. ( |3.35| ) with n = 4. This confirms that the cut- 
constructible terms of the amplitude do follow the naive factorisation of eq. ( 6.39| ) 

The rational terms of eqs. ( [5. 17] ) and ( |5.2| ), are each apparently singular in this limit. 
However, careful combination reveals the soft behaviour, 

i?4(</',l-,2+,3-,4+) + Ci?4(0,l-,2+,3-,4+)^C°^^(o)(l-,2+,3-,4+). (6.46) 

This is similar to the soft limit found in ref. [20, 64] for the MHV amplitudes with adjacent 
negative helicities, but, as anticipated in ref. [50], is not consistent with the naive limit of 
eq. ( p9| ). 
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7. Conclusions 



Previous analytic calculations of 0-amplitudes at one-loop with arbitrary numbers of gluons 
are the adjacent minus (/>-MHV [20], the all minus [51], and the finite all plus and single 
plus [50] i?!)-amplitudes. Higgs amplitudes produced by the effective interaction between 
Higgs and gluons induced by a heavy top quark loop, may be constructed from the sum of 
a i;^-amplitude and its parity conjugate cj)^ . In this paper, we have extended the calculation 
of one-loop MHV ^-amplitudes to include the general MHV configuration. 

One-loop amplitudes naturally divide into cut-containing, C„, and rational, parts. 
As in Ref . [20] , we used the double-cut unitarity approach of ref . [4] to apply the one-loop 
MHV rules to derive all the multiplicity results for the cut-constructible contribution Cn- 
In this paper we also used the spinor integration technique of ref. [6, 7] to determine Cn, 
finding complete agreement between the two methods. We found that the cut-constructible 
terms had a natural decomposition in terms of the pure glue MHV amplitude, we discovered 
that the new diagrams which arose as a result of the (j) interaction could be easily described 
by the basis functions used in the construction of the pure glue result. An explicit formula 
for the cut-constructible part of the (/>-amplitude are given in eq. ( |3.50|) , with the gluonic, 
fermionic and scalar contributions given in eqs. ( FsTl) , (|3^ ) and (|]5|). 

The rational terms have several sources - first the cut-completion term CRn which 
eliminates the unphysical poles present in Cn, second the direct on-shell recursion con- 
tribution R^, third the overlap term On and finally from the large z limit of the cut 
completion terms Inf CRn- Explicit formulae for each of these contributions are given in 
eqs. ( IS^STp , ( plj ), (|]2|) and respectively. 



The four gluon case is worked through in detail, and an explicit solution for the (j)- 
amplitude with split helicities, ^4^'*((/>, 1~, 2"^, 3~, 4"*"), together with instructions for how 
to assemble the Higgs amplitude A^p{H, 1^,2^,3^,4+) are given in section ||. Numerical 
results for this amplitude have previously been obtained in Ref. [61]. We have checked our 
analytic expressions in the limit where two of the gluons are collinear, in the limit where 
the (j) becomes soft and against previously known results for up to four gluons. 
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A. Evaluation of the Q, T and S functions. 
A.l + + 

The function Q is defined in eq. (|3.22| ). Using the Schouten identity, we can rewrite it as, 

Q(i,i^ + 1) = g{t,j) + g{t + 1, j + 1) - g{i + 1, j) - g{i,j + 1), (a.i) 

where the function g{i,j) is given by. 



Clearly, 

If i 7^ J then 

where 



g{i,i) = l. (A.3) 



N{P,i,j)= P^p, ■ pj - 2P ■ piP ■ pj . (A.5) 

g is now written in terms of scalar integrals so we can directly use the results of van 
Neerven [54] to perform the phase space integration: 

N{P,Pl,P2) 



J d^LIPS(-/i,/2,i^)^ 



2ism{7re)fs'^\P'r2Fi[l,-e;l-e;^^ -\ (A.6) 

(47r)^e^ V N{P,pi,p2)J 

d^UPS{-h,l2,P)^^^f^ = --^2isin(7re)/i2^|p2|- (A.7) 

I d^LIPS(-/i,/2,P) =-^^^^p^--^2isin(^e)^2e|p2|-. (^.8) 
where the factor cr is given by, 

,,„2 r(i + 6)r2(i-e) 

cr = (4vr) _ 2,) • (^-9) 

The final integration is over the z variable. However, the only dependence on z appears 
through the quantity Pi,n^ so it is convenient to make a change of variables, 

^ = JfflL (A.10) 

to produce a dispersion integral that will re-construct the parts of the cut-constructible 
amplitude proportional to {si^n)~^, 

'\2 



J ;^^2isin(^6)|Pr^ = 2ni{-P'r. (A.ll) 



^Through a suitable choice of 77, one can always ensure that N(P,pi,p2) is independent of z [4] 
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We define the function G{i,j) to be the reconstructed contribution after integration over 
phase space, and after performing the dispersion integration, 

G{i,j) = 1^1 d''UFS{-h,l2,P) G{i,j). (A.12) 

Exphcitly, we find that in the channel 

^ 5 (^)" (- (^-^ ' - « MB) ^ T^) 

The terms associated with triangle and bubble contributions will always cancel in the 
summation of the G{i,j) leaving only the contributions from the hyper geometric function 
as one would expect. 



A. 1.1 Spinorial Integration 

Let us show how the function G can be computed via spinorial integration. It is convenient 
to rearrange the integrand by applying different Schouten identities from the ones used 
above, so that 

gii,i + + 1) = -T{i,j) + r{i + 1, j) - r{i,j + 1) + r(i + 1, j + 1) (a.m) 

where 

r(„-).iiMl«^ ,A.i5) 

By using momentum conservation, I2 = P + /i, one can rewrite T{i,j) in terms of ^i, 

{ih) {j\P\h\ 
Then, one uses the rescaling in eq. (|3l0| ), so that, 

The above expression is the integrand of the double-cut integration, defined as, 

G'iiJ) = J dLIPS(4) r(i, j) . (A.18) 

By substituting the parametrization of dLIPS^^-* given in eq. ( |3.9| ), one has 



{£\p\£] J \ {e\p\e]) t {ie){j\p\£] 

{£ d£)[£ d£] P^ 



,2 (ij) 



{i£) {£\P\i] {j\P\£] 
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where the t-integration has been performed trivially. Before carrying through the spinor in- 
tegration, we introduce a Feynman parameter to combine the two denominators depending 
on 1^) 

G'iiJ) = [\x [ {i di)[i d£] P^ ^'f^j^'] (A.20) 

^ (27r)4 7 \ J {c\p\e] {i\R\ef 

where 

fi = x^i + {l-x)^ . (A.21) 
Integrating-by-parts in using the idenity, 

we obtain, 

The integration on \P\ can be performed by Cauchy's residues theorem, by taking the 
residues at the two poles, \i\ = f'\j) and \P\ = 

^^'^^ - {2^Y io I R-" {j\PR\j) {j\PR\j) 01 P RP\3] ] ^ ' 

Inserting the definition of ft in terms of x (paying attention to that is quadratic in x), 
we can perform the parametric integration, and by using some spinor identities, find that 

= / .LIPS(^) r,„-) = ^ in (1 - (^^^^fL_) (A.25) 



(2it)' V {2P ■ Pi)(2P ■ pj) j ■ 

which corresponds to the (discontinuity of) the double-cut of (the finite part of) the one- 
loop box function. 

A.2 ^(i,i + l,j,j + l) 

The function !f is defined in eq. ( |3.30| ). Again we define 

J'{i,i + 1, j, j + 1) = + l,i + 1) - Hi + 1, j) - Hhj + 1) (A.27) 

with, 

^ {^m) {jm){lh) jU,) 
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Then after using the Schouten Identity twice this can be written as, 



_ {il){im){jm){iei) , {im){jm){lj){l£2) , {i I) {i m) {j 1) {j m) {h £2 



{ij) (i^i)(lm)2 
Promoting to traces 



{ij) {j£2){lmY 



(A.29) 



tr_(l,i, j,m)tr_(l,£i,i,m) tr_(l, j, m) tr_(l, £2, j, H 



+ 



tr- (1, i, j, m) tr- (1, i, i, m) tr_ (j, i, I1J2) 
sj^slJ2i2.pj)i2£i.p.) 



(A.30) 



Which we recognise as two hnear triangles and a box function similar to those in ^. If we 
commute £2 and j in the final term we can get something which looks like eq. (A. 4). 

tv_{j,i,£i,£2) , tT-{j,£i,i,£2) 



{2£i.pi){2l2.Pj) 



1 



{2£i.pi){2l2.pj) 



(A.31) 



The first term will cancel the bubbles which arise in the calculation and the remaining 
terms are triangles and boxes. However, since the coefficients of J- depend on i and j there 
will no longer be a cancellation between the four terms. This is important in controlling 
the IR divergences of the amplitude, the triangle pieces are needed to cancel off the IR 
poles coming from the box functions. After performing Passarino-Veltman reduction on 
the tensor integrals and performing the dispersion integrals we find. 



p2 



tr_(l,i, j, m) tr_(l,j,i,m) 



„2 2 



1 



Fi l,-e;l-e; 



. Pi ■PjP'^ 



N{P,i,j) 



1 



/ tr_(l,i, j,m)tr_(l,P,i,m) 1 tr^(l,j,i,m) tr_(l, P, j, m; 

V Si.sj^ 2(P.k) 2iP.p,)Jl-2e 



(A.32) 



A. 2.1 Spinorial Integration 

Alternatively the function can be computed via spinorial integration. Using momentum 
conservation, I2 = P+h, and the rescaling in eq. ( 3.10 ), one can rewrite of eq. ( |A.2^ ) 

in terms of £, and t, 

^ {rm) ijm) {£!) {l\P\£] 
{i£){mlf{j\P\£] 

The above expression, which turns out to be independent of t, is the integrand of the 
double-cut integration, defined as, 

F'iiJ) = J dLIPS(4) J^iiJ) . (A.34) 
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By substituting the parametrization of dLIPS*-^-* given in eq. ( |3.9D , and performing the 
phase-space integration with spinor- variables, one finds, 

j,,^^^^^ = _ H^Jl\rn]HJ^l\m] ^,^^^^^ ^ 

^ij ^Im 

27Ti f {m\j il\m]{m\i Pl\m] \ 

where G'{i,j) was given in eq. ( |A.26D . We remark that the term proportional to G'{i,j) 
corresponds to the (discontinuity of) the double-cut of (the finite part of) the one-loop 
box function; while the rational part of eq. ( |A.35| ) corresponds to the discontinuity of 
logarithmic functions associated with a combination of 2-point and (Im- and 2m-) 3-point 
functions. 

A.3 5(i,i + l,j,j + l) 

The final pieces of the amplitude, associated with the propagation of scalar particles around 
the loop, are the most complicated. The function S is defined in eq. ( |3.31 ). In a similar 
fashion to the gluonic and fermionic pieces we define, 

S{i,i + + 1) = S{i,j) + S{i + 1, i + 1) - S{i + S{i,j + 1) (A.36) 

with 

.^ {lhf{U2f{mh){m£2) {im) {j m) 

on, J) = 2 K (A. 37) 

{lm)Uhi2f{iii){ji2) 
After using the Schouten Identity the above can be reduced to a scalar box and third rank 
triangles which can be solved via Passarino-Veltman reduction generating, 

^ '"tr_(l, i, j, m)^ tr_(l, j, i, m)^ ^ ^ ^ ,v,.d2 



2 I _p2 



4 4 

f 1 /tr_(l, i, j, m) tr_(l, P, i, m)^ 1 

1 /tr_(l, i, j, m)^ tr_(l, P, i, m)^ 1 

tr_(l,i,j,m)2tr_(l,P,i,m)tr_(l, j,i,m) 1 _^ ■ ^ 



44m (2K-^) 7Jl-2eJ 

(A.38) 



A. 3.1 Spinorial Integration 



Alternatively the function S can be computed via spinorial integration. By using momen- 
tum conservation, I2 = P + h, and the rescaling in eq. ( 3.101 ), one can rewrite eq. ( |A.37| ) 
in terms of £ and t, 

^ ^2 (im) (jm) (m£) {llf {m\P\e\{l\P\lf 
P^{ii) {ml)^{j\P\£] 
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By substituting the parametrization of dLIPS*-^-* given in eq. ( |3.9D , and performing the 
phase-space integration with spinor- variables, one obtains, 

S'iiJ) = j dLIPS(") 

(m| i i 1 |m]^ (m| j i 1 |?ti]^ , 27rz J (m| 1 z j |m] (m| i P 1 |m] 

^ (m| i P 1 ^ (m| 1 i j |m] (m| i P 1 \m] 



3 (i|P|i]^ 2 Sij {i\P\i] 

(ml 1 i 7 Iml (ml i i 1 Iml \ , , 1 / . . s 

+ — [2 \ , (A.40) 



where G'{i,j) was given in eq. ( |A.26D . We remark that the term proportional to G'{i,j) 
corresponds to the (discontinuity of) the double-cut of (the finite part of) the one-loop 
box function; while the rational part of eq. ( |A.4C| ) corresponds to the discontinuity of 
logarithmic functions associated with a combination of 2-point and (Im- and 2m-) 3-point 
functions. 



B. Scalar integrals 

The one-loop functions that appear in the all-orders cut-constructible contribution Cn given 
in section 3 are defined by, 



Fl^'{P\Q^-s,t) = ^ 



2 \ <^ 



p2 

2 \ <^ 
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and 



Bub(s) 
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e(l - 2e) V-s 



s ^ 

2 \ e 



(B.4) 
(B.5) 
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